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Abstract 

We consider the evolution of an interface generated between two immiscible incom- 
pressible and irrotational fluids. Specifically we study the Muskat and water wave prob- 
lems. We show that starting with a family of initial data given by (a, /o (<*)), the interface 
reaches a regime in finite time in which is no longer a graph. Therefore there exists a 
time t* where the solution of the free boundary problem parameterized as (a, f(a,t)) 
blows-up: ||9 a /||ioo(t*) = oo. In particular, for the Muskat problem, this result allows 
us to reach an unstable regime, for which the Ray leigh- Taylor condition changes sign and 
the solution breaks down. 

1 Introduction 

Here we study two problems of Fluids Mechanics concerning the evolution of two incompress- 
ible fluids of different characteristics in 2D. We consider that both fluids are immiscible and 
of different constant densities p 1 and p 2 , modeling the dynamics of an interface that separates 
the domains (t) and f2 2 (i). That is, the liquid density p = p(x,t), (x,t) € ]R 2 x M + , is 
defined by 

o(x t)-l p1 ' xeQl ^ (i) 

p^-\ p 2 , xen 2 (t) = R 2 -n 1 (t), [i) 

and satisfies the conservation of mass equation 

* + t,-V, = 0, 
V-v = 0, 

where v = (vi(x,t),V2{x,t)) is the velocity field. With a free boundary parameterized by 

dQP(t) = {z(a,t) = (zi(a,t),z 2 (a,t)) : a € R}, 
we consider open curves vanishing at infinity 

lim (z(a,t) — (a,0)) = 0, 

a— >oo 

or periodic in the space variable 

z(a + 2kir,t) = z(a,t) + 2/^(1,0). 
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The scalar vorticity, V • v, has the form 

V" 1 • v(x,t) = uj(a,t)5{x - z(a,t)), (3) 
i.e. the vorticity is a Dirac measure on z defined by 

<V ± -v,r/>= / u(a, t)r](z(a, t))da, 
Jr 

with r](x) a test function. The system is closed by using one of the following fundamental 
fluid motion equations: 

Darcy's law 

Hv = -Vp-gp(Q,l), (4) 

K 

or 

Euler equations 

p(v t + v ■ Vv) = -Vp-gp(0,l). (5) 

Here p = p(x, t) is pressure, g gravity, p viscosity and k permeability of the isotropic medium. 

The Muskat problem [18] is given by (jl!2l4p which considers the dynamics of two incom- 
pressible fluids of different densities throughout porous media and Hele-Shaw cells [221 H5] • In 
this last setting the fluid is trapped between two fixed parallel plates, that are close enough 
together, so that the fluid essentially only moves in two directions |15j . 

Taking p 1 = 0, equations (|1I2I3I5P are known as the water waves problem (see [2] and 
references therein), modeling the dynamics of the contour between an inviscid fluid with 
density p 2 and vacuum (or air) under the influence of gravity. 

Condition ([3]) (deduced by (j4]), assumed for (j5])) allows us to write the evolution equation 
in terms of the free boundary as follows. One could recover the velocity field from ([3]) by 
means of Biot-Savart law 

v(x,t) = V L A~\V L -v)(x,t) = ±- [ (*~ u(a,t)da, 

2vr J R \x - z(a,t)\ z 

applying the delta measure with amplitude to. Taking limits on the above equation approach- 
ing the boundary in the normal direction inside fi- 7 , the velocity is shown to be discontinuous 
in the tangential direction, but continuous in the normal, and given by the Birkhoff-Rott 
integral of the amplitude u along the interface curve: 

BR(z,u)(a,t) = ±PV [ ^ -^^ (MdP- 
2vr J R \z(a,t) - z{/3,t)\ z 

Above PV denotes principal value. It yields the curve velocity from which one can subtract 
any term c in the tangential without modifying the geometry of the interface 

zt(a, t) = BR(z, u)(a, t) + c(a, t)d a z(a, t). (6) 
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Understanding the problem as weak solutions of (|1I2|4|) or (jl!2|3l5p . the continuity of the 
pressure on the free boundary follows. Therefore, taking limits in Darcy's law from both 
sides and subtracting the results in the tangential direction, it is easy to close the system for 
Muskat (in this paper we consider two fluids with the same viscosity): 

w(a,t) = -{p 2 - P l )—d a z 2 {a,t). (7) 
A 1 

In a similar way for water waves, Euler equations yield 

uj t (a,t) = -2d t BR(z,u)(a,t) ■ d a z(a,t) - d a (-^ — rj )(«,*) + d a (cLu)(a,t) 
+ 2c(a, t)d a BR(z, u>)(a, t) ■ d a z(a, t) + 2gd a Z2(a, t). 



Then, the two contour equations are set by (|6|7p and ([6f 

For these models the well-posedness turns out to be false for some settings. Rayleigh 
|21j and Saffman- Taylor [22] gave a condition that must be satisfied for the linearized model 
in order to exist a solution locally in time: the normal component of the pressure gradient 
jump at the interface has to have a distinguished sign. This is known as the Rayleigh- Taylor 
condition. It reads 

a(a,t) = -(Vp 2 (z(a,t),t) - Vp 1 (z(a,t),t)) ■ d^z(a,t) > 0, 

where Vp>(z(a, t), t) denotes the limit gradient of the pressure obtained approaching the 
boundary in the normal direction inside £l 3 (t). 

An easy linearization around a flat contour (a, f(a,t)), allows us to find 

ft = \h{u) 

where H is the Hilbert transform which symbol on the Fourier side is given by H = — i sign (£). 
The equations 

uj = —(p 2 — p 1 )—d a f, (linear Muskat) 

A* 

uit = 2gd a f, (linear water waves) 

show the parabolicity of the Muskat problem when the denser fluid is below (p 2 > p l ) and 
the dispersive character of water waves. 

1. There is a wide literature on the Muskat problem and the dynamics of two fluids in 
a Hele-Shaw cell. There are works considering the case of a viscosity jump neglecting the 
effect of gravity |23| 128] . Local-existence in a more general situation (with discontinuous 
viscosity and density) is shown in [7] and also treated in [I]. A different approach to prove 
local-existence can be found in |10| for the setting we are considering in this paper. The 
Rayleigh- Taylor stability depends upon the sign of (p 2 — p l )d a zi(a,t) indicating that the 
heavier fluid has to be below in the stable case. If the lighter fluid is below, the problem 
have been shown to be ill-posed jlUj . Global-existence results for small initial data can be 
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found in 0, [2U [Ml QUI E] . For large initial curves and parameterized by (a, f(a,t)), there 
are maximum principles for the L°° and I? norms of /, and decay rates, together with 
global-existence for Lipschitz curves if ||3 Q ,/||i / cx>(0) < 1 [TT| [5]. 

2. The water waves problem have been extensively considered (see [H [2] and references 
therein). For sufficiently smooth free boundary, the Rayleigh- Taylor condition remains pos- 
itive with no bottom considerations [25] . a fact that was used to prove local-existence [25]. 
The Rayleigh- Taylor stability can play a different role for the case of non "almost" flat bot- 
tom [T7]. Recently, for small initial data, exponential time of existence has been proven 
in two dimensions [26], and global-existence in the three dimensional case (two dimensional 
interface) OH [27]. 

The results that we announce below will be given in the forthcoming papers [3] and [4]. 



2 Rayleigh- Taylor breakdown for Muskat 

This section is devoted to show the main ingredients to prove the Theorem 12.11 below. We 
consider the function 

F(zo)(a,P)= ^ ™ (9) 

\z {a) - z (a - P)\ 2 

if F{zq) G L°°(IR 2 ) then the curve zq satisfies the arc-chord condition. We say that the 
Rayleigh- Taylor of the solution of the Muskat problem breaks down in finite time if for initial 
data zq satisfying a(a,0) = (p 2 — p 1 )d a zi(a,0) > there exists a time t* > for which 
a(a,t*) is strictly negative in a nonempty open interval. 

Theorem 2.1 There exists a non-empty open set of initial data in H 4 , satisfying Rayleigh- 
Taylor and arc-chord conditions, for which the Rayleigh- Taylor condition of the solution of 
the Muskat problem (El IB breaks down in finite time. 

After choosing the appropriate tangential term and a integration by parts, the contour 
equation reads 

P 2 -P X py f (zi(a,t) - Zl ((3,t) ) 
2tt J \z(a,t) - z((3,t)f 



zt(a,t) = P -^PV I ^0-(d a z(a,t)-d a z(M)dP. (10) 



The steps of the proof: 

1. First, for any initial curve zo(a) = z(a, 0) in H A that satisfy R-T 

{p 2 - p l )d aZl {a^) >0 

and the arc-chord condition then the solution to the Muskat problem z(a, t) becomes analytic 
for < t < T. Moreover, z(a, t) is real analytic in a strip 

S(t) = {a + i(: |C| < ct} 

for t € (0, T) where c depends only on 

d a zi(a,0) 



inf (0) = inf 



\d a z(a,0)\ 2 - 
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The proof follows by controlling the quantities extended on S(t): 

f ( . ., ,x d a zi(a + iC,t) 

f{a+lC > t) = \d aZ (a +l c,t)r 

F(z)(a + i(,(3,t) by using ([9]) and norms 

11-^0) ||z,°° (s) (*) = SUP \F(z)(a + i(, 

a+i(eS{t),l3£j 

IMIl 2 (S)W = / \ z i a ^ i c tit)\ 2 da, 



|2 



(t) = \\4h iS )(t) + ^2 J^Kz(a ± ict,t)\ 2 da, 



for j € 



inf(t) = inf »(/)(a + »C). 

a+iCG5(t) 



Then the quantity 

IWI*r(*) = 11*11*4(5) W + ||F(z)|| Loo(s) (t) + l/(inf(t) - c - K\\Q(f)\\ H 2 {S) (t)) 

satisfies 

j t \\z\\ RT (t)<C\\z\\ k RT (t), 
for C, -fT and k universal constants. It yields 

\\4rt(o) 



\z\\RT{t) < 



(l-C\\z\\ RT (0)t) l/k ' 



providing control of the analyticity and T = l/(C||z||^ T (0)). 

2. Second, there is a lower bound on the strip of analyticity, which does not collapse to 
the real axis as long as the Rayleigh- Taylor is greater than or equal to 0. Then there is a 
time T and a solution of the Muskat problem z(a,t) defined for < t < T that continues 
analytically into a complex strip if (p 2 — p l )d a z\ > 0, where T is either a small constant or 
it is the first time a vertical tangent appears, whichever occurs first. We redefine the strip 

S{t) = {a + i(: |C| < h(t),0 < h(0)}, 

and the quantity ||^||| = IMI^Wg) + II^ 71 ( z )IIl oo (5) with this new S(t). For a h(t) decreasing 
(the expression of h(t) is chosen later), we consider the evolution of the most singular quantity 



± 

Taking a derivative in t one finds 
d 



/ \d*z(a±ih(t),t)\ 2 da. 



dt ^2j Kz(a±ih(t))\ 2 da<^-Yl J Md*z)(a±ih(t))-d*z(a±ih(t))da 
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-10h'(t) J A(d^z) (a) ■d*z(a)da + 2^3? J d*z t (a ± ih(t)) ■ d*z(a ± ih(t))da. 
Estimating in a wise way one obtains 

A£/|^ ±! M t »^ < owiw - io A '( t )/A(^,)( ).^ 

+(C||«|||(t)/i(t) + — h'(t)) J k(diz){a±ih{t))-d&{a±ih{t))da. 
Therefore, choosing 

fe(t) = /i(0)exp(-10C f \\z\\ k s (r)dr) 
J o 

eliminates the most dangerous term. The other terms are easily controlled, giving finally 

|E/ \dUa±ih(t))\ 2 da<C\\z\\ k s (t). 

This allows us to reach a regime for which the boundary z develops a vertical tangent at time 
T. 

3. Third, it is shown the existence of a large class of analytic curves for which there exist 
a point where the tangent vector is vertical and the velocity indicates that the curve is going 
to turn up and reach the unstable regime. That is 

a. d a z\{a) > if a ^ 0, b. d a zi(0) = 0, 

c. d a z 2 {0)>0, d. 9 a «i(0) < 0, 

for analytic functions z\{a) and z 2 {a) such that z(a) satisfies the arc-chord condition. Here 
we consider an open curve vanishing at infinity (being analogous in the periodic case). We 
assume that z(a) is a smooth odd curve satisfying the properties a, b and c. Differentiating 
the expression (llOp for the horizontal component of the velocity it is easy to obtain 

,q Vl \r a \ = f°° (gafl(g) ~ d a zi{a- ft)) 2 + (zi(q) - zi(q-/3))(^zi(q) - dlz^a-p)) ^ 
1 J-oo \z{a) - z(a-/3))\ 2 

(zi{a)-zi{a-f3))(d a zi{a)-d a zi(a-p)) _ ^ a _^| 4 rf /3. 

At a = it yields 



i) - z 2 (f3)(d a z 2 {<d) - d a z 2 i 

Integration by parts provides 



\m\ 2 1 Loo k(/?)i 2 dp 



-oo 

Therefore it is easy to obtain that 



(d aVl )(0) = 4d a z 2 (0) I TJl^A ' d a z x {P)dp. (11) 







zi{P)z 2 {P) 

W)\ 

Expression (jlip allows us to determine the sign of (d a i>i)(0). One could take 

P 3 



zi{P) 



(1 + /3 2 ) 



and construct the function z 2 {P) in the following way: Let ft, ft and ft be real increasing 
numbers. We pick z 2 (P) < c < for ft < P < oo and z 2 (P) a smooth function with the 
following properties 

a. z* 2 {P) is odd, b. (dpx%)(0) > 0, 

c. z* 2 (p) > if P E (0,/3i), d. z* 2 (p) < if P E (ft, 
For z 2 (P) = bz 2 (P), < P < ft and 6 > 0, the velocity satisfies 

(d a vi){0) < 4{d a Z2){0) \ z (p)\4 d a z nP) d P + \z(P)\ 4 daZl ^P> d P 

= 4(8 ^ )(0) (T (WTiiw 8 " 2l( ' 3) * + A ) ■ 

where A < 0. The constant 6 large enough yields (d a vi)(0) < 0. Rectifying the curve on the 
interval [ft, ft] it is easy to obtain a smooth curve. Finally, convolving with the heat kernel 
the vertical component, the curve z(a) is approximated by an analytic one. 

4. Fourth, with the initial data found in 3. and no assumption on the R-T condition, 
we use a modification of Cauchy-Kowalewski theorems |19^ [20] to show that there exists 
an analytic solution for the Muskat problem in some interval [— T, T] for a small enough 
T > 0. Here we are forced to change substantially the method in [24] because in this case the 
curve can not be parameterized as a graph, so we have to deal with the arc-chord condition. 
Then, with {X r } r> Q a scale of Banach spaces given by real functions that can be extended 
analytically on the complex strip S r = {a + iC, € C : \Q\ < r} with norm 



± 



and z°(a) a curve satisfying the arc-chord condition and z° (a) E X ro for some tq > 0, we 
prove the existence of a time T > and < r < ro so that there is a unique solution to the 
Muskat problem in C([—T,T];X r ). This allows us to find solutions that do not satisfy the 
R-T but shrink the strip of analyticity. We extend equation (IIOD as follows: 

z t (a + i(, t) = G(z(a + i(, t)), 
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with 



G(z)(a + <, *) = y R^cF^T<^F (aaZ(a + 1 ° " 9aZ(a + * c ~ 

For < r' < r and the open set O in <SV given by 

O = {z E X r : \\z\\ r < R, \\F(z)\\ Loc(Sr) < i? 2 }, (12) 

the function G for G : O — >• X r / is a continuous mapping and there is a constant Cr (depend- 
ing on only) such that 

\\G(z)\\ r , <^||4, (13) 

||G(z 2 ) - g^ 1 )^ < ^lk 2 - (14) 

and 

sup |G(z)(Q + iC)-G'(z)(a + iC-/3)| <Cr|/3|, (15) 

for z, z^ E O. For initial data 2° G X ro satisfying arc-chord, we can find a < Tq < ro and a 
constant i?o such that ||z°|| r ^ < -Ro an d 

(z (a + - C) _ z {a + K _ p)) 2 + (z (a + - C) _ z (a + K _ p)) 2 ^ 



3 2 



Rl 



for a + iQ E 5 r / . We take < r < 7q and Rq < R to define the open set O as in (fT2j) . 
Therefore we can use the classical method of successive approximations: 



z n +\t) = z°+ f G{z n {s))ds, 
Jo 

for G : O — > X r i and < r' < r. We assume by induction that 



"{t)<R, and \\F(z k )\\ Loo(Sr )(t) <R 



z 



for k < n and < t < T with T = min(T J 4, Tqk) and Tex the time obtaining in the proofs 
in [19] and [20]. We get ||z n+1 || r (t) < R that follows using (|13|14p . The time T A is to yield 
||F(2 n+ )||x,oo(5 r ) (t) < R. Then, using the induction hypothesis and ([To]) we can control the 
quantity taking < T A < (Rq 2 - R- 2 ){C 2 R + 2 j R Cr)~ 1 - 

5. Fifth, all the results above allow us to prove that there is a non empty set of initial 
data in H A satisfying the arc-chord and R-T conditions such that the solution of the Muskat 
problem reaches the unstable regime: the R-T becomes strictly negative on a non-empty 
interval. We pick initial data as in 3. We apply the local-existence result in 4 to get an 
analytic solution z{a, t) on [— T, T]. Then we consider a time < 5 < T and a curve wf (a, t), 
solving the Muskat problem with initial datum z(a, —5) + rrj(a). The function rjg has a small 
H A norm, i.e. 

\\UI(;-S) - Z(;-S)\\ m = < 8. 
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The time 5 is small enough so that w|(a, —S) satisfies R-T: (p 2 — p 1 )d a (uj$)i(a>, —5) > 0. Then 
we apply the local existence result in 1. that wf (a, t) becomes analytic for some time —6 < t. 
With 2. we assure the existence and analyticity of the solution even if d a {u)f)\{a,t) < for 
some time t. Then, we show that both solutions are close in the H topology as time evolves. 
We can apply to w| the local-existence result in 4. if it is needed. Then, with 5 and e small 
enough we find the desired result. 

3 Turning water waves 

In this section we prove for the water wave problem (p 1 = and (|1I2I3I5P ) that with initial 
data given by a graph (a,/o(a)), the interface reaches a regime in finite time where it only 
can be parameterized as z(a, t) = (z±(a, t), z%[a.^ t)), for a£l, with d a zi(a, t) < for a € /, 
a non-empty interval. Therefore there exists a time t* where the solution of the free boundary 
problem reparameterized by (a, /(a, i)) satisfies ||/ q ||l°°(**) = °o- 

Theorem 3.1 There exists a non-empty open set of initial data (a,/o(a)), with /o £ H 5 , 
such that in finite time t* the solution of the water waves problem (p 1 = and CTIM ) 
given by (a, f(a,t)) satisfies \\fa\\L°° (t*) = oo. The solution can be continued for t > t* as 
z(a,t) with d a zi(a,t) < for a £ I, a non-empty interval. 

In order to prove this theorem we consider a curve z*(a) € H 5 with the same properties 
as in point 3. of previous section. Then, we pick z(a,t*) = z*(a) and u(a,t*) = d a z^(a) 
as a datum for the initial value problem. It is easy to find the same properties for the 
velocity, since the tangential direction does not affect the evolution. Picking the appropriate 
c(a,t) and applying the local-existence result in [8], there exists a solution of the water 
waves problem with z(a,t) G C([t* - S,t* + 5};H 5 ), u(a,t) G C([t* - S,t* + 5};H*) and 
5 > small enough. Then, the initial datum (zo(a),uJo(a)) = (a, fo(a), uo(a)) is given by 
(z(a,t* - 5),u;(a,t* - 5)). 

4 Muskat breakdown 

In this section we show that there exists a smooth initial data in the stable regime for the 
Muskat problem such that the solution turns to the unstable regime and later it breaks down. 
The outline of the proof is to construct a curve in the unstable regime which is analytic except 
in a single point. We show that as we evolve backwards in time the curve becomes analytic 
and is as close as we desired (in the H k topology with k large enough) to the curve from part 
3. of section 2. 

Here we will work in the periodic setting and will consider the equation 



MC,*)= / 



sin(>i(C,t) - Z!(w,t)) 



(d^z((, t) — d(z(w, t))dw, 
(17) 



J w 



er+(t) cosh(z 2 (C,t) - z 2 {w,t)) - cos(^(C,t) - zi («>,*)) 



where ( £ fi(t) 



Q(t) = {(G C/2kir : |3?C| < h(Jftz,t)} 
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h(x, t) is a positive periodic function with period 2tt and smooth for fixed time t and 

r ± (t) = {C G C/2/cvr : ( = x + ih(x, t)}. 

This equation is equivalent to (jl!2l4p for holomorphic functions. 
In order to prove the result we will need the following theorem: 

Theorem 4.1 Let h(x, t) be a positive, smooth and periodic function with period 2ir for fixed 
time t G [to — 6, to]. Let z(x,to) be a curve satisfying the following properties: 

• z\(x,to) — x and Z2(x,to) are periodic with period 2tt. 

• z((,to) is real for £ real. 

• z((,to) is analytic in C G Q(to)- 

• z{C,,t) G H k (T±{to)) with k a large enough integer. 

• Complex Arc- Chord condition. 

| cosh(z 2 (C,t ) - Z2(w,t )) - cos(^(C,to) - *i(w,to))| > [\MC -w)\\ + |9f(C - w)|] 2 , 
/or (, w£ S7(io)j where \ \x\ \ = distance(x,2kir). 

• Generalized Rayleigh- Taylor condition: RT(^,to) > 0, where 

RT{C " = s ( (^(c^+( ( ^(c.o)' (1 + * ft(SC ' 

+s({p-V.[ r, 77 , sin( ; l(C '' ) ; ^ A» + .9,MC, 1 

\[ 7^r+(i) cosh(2 2 (C,t) - z 2 {w,t)) - cos(zi(C,t) - 2i(w,i)J J 

x(l + ia :c / l (KC,i)) _1 ) 

Then, for small enough 5, there exist a solution for the equation JJ7| j in the time interval 
t G [to — 5, to] satisfying 

• z\(x,t) — x and Z2(x,t) are periodic with period 2ir. 

• z(C, t) is real for £ real. 

• z(<^,t) is analytic in £ G Sl(to) • 

• z((,t) G .ff fc (r-|-(t)) wit/i A; a large enough integer. 

Now, let z(x,t) be the solution of the Muskat problem with z(x, 0) = where z®(x) 

is the particular initial data from part 3. of the section 2. We shall define this solution as 
the unperturbed solution. Let us denote the Rayleigh- Taylor function 

11 ' ; ~ (^iM)) 2 + fe 2 M)) 2 ' 

Notice the minus sign in the right-hand side of the previous expression. One can check the 
following properties of this Rayleigh- Taylor function: 
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1. (Ti(-,t) is analytic on {x + iy : iGT, \y\ < q,} with |<rj(a; + iy,t)\ < C, for all x + iy 
as above and for all t < [0, r]. 

2. cj?(0, 0) is real for x G T, t G [0, r]. 

3. cjj 1 has a priori bounded C k ° norm as a function of (x,t) G T x [0, r] (fco large enough). 

4. of (0,0) = 0. 

5. 5 x .cj0(0,0) = 0. 

6. <9 2 ct?(0,0) = -02 < 0. 

7. <Vi(0,0) = a > 0. 

In this setting we define the following weight functions 

fc(sc,t) = ^- 1 (r 2 -t 2 ) + (^- 1 -(r-t))sin 2 (|) fortG[r 2 ,r]. (18) 

h{x,t) = - (A- 1 r 2 + A" 1 sin(|)) + A~ 2 Tt + At sin (|) tG[0,r 2 ]., (19) 

with x G T. First we choose the parameters A large enough and then r small enough, then 
one can show that 

<7?(x,t) +d t h(x,t) - A^h{x,t) > ct 2 for x G T, t G [r 2 ,r] (20) 

and 

a?(x,t) +d*£(x,i) -^H(x,t) > -A _2 r for x G T, t G [0,r 2 ]. (21) 

The inequalities (|20|) and (|21|) are one of the main ingredients of the proof of the following 
results 

Theorem 4.2 Let z(x,t) be a solution of the Muskat equation in the interval t G [0, r]. Let 
h(x,t) and h~(x,t) as in the expressions U8\) and U9\) and k a large enough integer. Assume 
that z(x, t) satisfies 

• z\(x,t) — x and Z2(x,t) are periodic with period 2ir. 

• z(C)i) is real for £ real. 

• z(C)i) *s analytic in £ G f2(i). 

• ^(C)i) £ -f^ fc (r±(i)) wfi fc a Zan/e enough integer. 

• Complex Arc-Chord condition. 

| cosh(z 2 (C,t) - z 2 K*)) - cos(z! (C,t) - 2i > [||»(C - w)|| + |9f(C - w)\} 2 , 
for C, w G H(t). 
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Here in the definition of £l(t) and T±(t) we use h(x,t) if t G [t 2 ,t] and h~(x,t) if t G [0,r 2 ]. 
Then 

2 



Id 
2~dt 



u>er + (t) 



d£z(C, t) - d£z((, t) d$l( > -C{A)\\ 



iftG[r 2 ,r] 



and A < r 50 . 
In addition 



j d k c z((,t)-d k c z((,t) 

Jwer + (t) 



d$t( < A 2 



ld_ 

2di 



d%z((, t) - d%z((, t) dStC > -C{A)t- 1 X\ 



lft£[0,T 2 



L 



w&r + (t) 



d k c z((,t)-d k c z((,t) 



d$t( < y 



and A < r 50 . 

This theorem implies that for all 7 > there is e > such that 



L 



w&T + (t) 



for t € [0, r] if 



d k ( z(C,T)-d k ( z((,T) 



MC < 7 



d$t( < e 



iwer+{t) 

and z(x,t) satisfies the requirements of the theorem. 

Lemma 4.3 Let z(x,t) be a solution of the Muskat problem satisfying the requirements of 
theorem {^..2) and close enough to the unperturbed solution int G [0, r]. Let h(x,t) and K(x,t) 
be as in U8\) and U9\) with a suitable choice of A and r. Then z(x,t) satisfies the generalized 
Rayleigh- Taylor condition in t G [0, r]. In particular the unperturbed solution satisfies the 
generalized Rayleigh- Taylor condition in t G [0, r]. 

Theorems (|4.ip . (|4.2f) and lemma (14. 3D allow us to achieve the desired result. Indeed we can 
choose a curve z(x,r) such that 



Cer± 



d k z(C,r)-d k z(C,r) 



d^C < e, 



with < e < eo ( £ o small enough), satisfying the generalized Rayleigh- Taylor condition by 
lemma (|4.3j) and satisfying the rest of the hypothesis of theorem (|4.ip . Since /i(0, r) = 0, 
z(x,t) is allow to be no analytic at x = (maybe z(»,r) G H k (T) but z(x,r) // fe+1 (T)) . 
By theorem (|4.ip there is a solution analytic in Q(t), for some interval t G [r — <5, r] 
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with small enough 5 and for all e. By theorem (|4,2p . we can choose e small enough in such 
a way that, by lemma (|4.3p . z(x,t — (5) satisfies the generalized Rayleigh- Taylor condition. 
Then we can go further the time r — 5. Iterating this argument, we find we can extend z(x,t) 
to be a solution of the Muskat problem, analytic in Q(t) for all t S [0, r] and as close as we 
want to the unperturbed solution. 
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